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We present numerical results and their analyses of finite-size effects of hadron masses for both quenched and 
full QCD calculations. We show that they are much larger for full QCD due to dynamical sea quarks and the 
associated breaking of Z(3) symmetry. We also argue that finite-size effects are non-negligible even for the largest 
lattice size simulation currently being made for a very small quark mass. 



1. Introduction 

In this report we summarize a comparative 
study of finite size effects for hadron masses in 
quenched and full QCD with the Kogut-Susskind 
quark action[l]. Our motivation for the study is, 
first, to examine whether the impression obtained 
from previous quenched studies that finite-size ef- 
fects are perhaps smaller for quenched QCD [2] 
is actually valid, and, second, to understand the 
origin of the difference if it exists, including the 
possible role of dynamical sea quarks in the large 
finite-size effects observed for full QCD[3]. 

Our two-flavor full QCD study[3] was carried 
out at (3 = 5.7 for the spatial lattice size in 
the range La = 0.7 - 1.8fm (L = 8 - 20) with 
a = 0.089(3)fm fixed by the p meson mass. For 
a comparative quenched study we chose (3 = 6.0 
where our previous calculation on a 24 3 x 40 lat- 
tice gave a = 0.105(3)fm[4], and made new runs 
for the sizes L = 6 - 16(ia = 0.63 - 1.7fm). We 
also added a full QCD run on an 8 3 x 16 lattice 
using quark boundary conditions different from 
those of Ref. [3]. 

2. Comparison of lattice-size dependence 

Our full and quenched QCD results for ir, p 
and N masses are compared in Fig. 1 where the 
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periodic boundary condition is imposed on sea 
and valence quarks in the spatial directions. It 
is evident in Fig. 1 that the magnitude of finite- 
size effects is much smaller for quenched QCD 
than for full QCD, especially below La ~ 1 fm. 
The size dependence is also significantly weaker 
for quenched QCD; assuming a power law b^m oc 
L~ a we find a « 1 — 2 for quenched QCD as 
compared to a « 2 — 3 observed for full QCD[5]. 

3. Origin of difference 

The difference observed in Fig. 1 can be un- 
derstood by a simple argument based on an ex- 
pansion in inverse powers of quark masses. The 
essential points of the argument are present in 
the literature, dating back to quite early times 
for the quenched case[7] and more recently for 
full QCD[5]. 

For meson propagators application of the ex- 
pansion for valence quarks leads to an expression 
of form 

c c 

where m va i is the valence quark mass and {■) de- 
notes the gluon field average (including the quark 
determinant in full QCD) . The second term comes 
from valence quark loops C wrapping around the 
lattice in the spatial directions (Polyakov-type), 
and the first from trivial ones (Wilson-type), with 
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Figure 1. Comparison of quenched[l,4] and full 
QCD [3,6] hadron mass data. Solid lines are fits 
of form m = + c/L 3 to full QCD data. 
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Figure 2. 7r and p masses for various quark 
boundary conditions at m q a = 0.01. 
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P(C) and W(C) the associated gauge link fac- 
tors. The sign factor a va i for the second term 
represents the spatial boundary condition taken 
for valence quarks; a va i = +1 or — 1 for the peri- 
odic (P) and the anti-periodic (AP) cases. 

It is obvious from (1) that hadron masses would 
be significantly shifted if (P(C)) takes a large 
value. This in fact is what happens in full QCD 
for small lattice sizes where breaking of center 
Z(Z) symmetry induced by dynamical sea quarks 
aligns Polyakov-type loops in some particular di- 
rection. In contrast, a sufficiently large ensem- 
ble of quenched configurations is invariant under 
Z(3), and hence (P(C)) = even for small lat- 
tices unless the winding number of the loop is 
an integer multiple of 3. Thus quenched hadron 
masses should suffer less from finite-size effects, 
as is indeed observed in Fig. 1. 

This line of argument can be pushed a step fur- 
ther and predict the sign of the mass shift for full 
QCD. If one imposes the anti-periodic boundary 
condition for sea quarks the Polyakov-type loops 
take positive values on the average, while for the 
periodic boundary condition they take negative 
values. Combined with the sign factor a va i due 
to valence quarks this means that the Polyakov- 
type terms cancel the Wilson-type contribution in 
(1) for the boundary conditions (P seC nPvai) and 
(AP aea , AP va i). Meson propagators will fall off 
faster for these cases, which would lead to an in- 
crease of meson masses. Similarly a decrease is 



expected for the opposite case of the (P sea , AP va i) 
and (AP sea , P va i) boundary conditions since the 
two terms add up. 

In Fig. 2 we plot meson mass data taken on an 
8 3 X 16 lattice with the four boundary conditions, 
which confirm these expectations. In the figure 
we also plotted quenched results for the P va i and 
AP va i boundary conditions on a 6 3 X 40 lattice. 
A small difference seen between the two cases for 
quenched QCD, negligible compared to those for 
full QCD, provides confirmation that Polyakov- 
type terms are almost cancelled out in quenched 
QCD as argued above. 

An interesting way to further check the cancel- 
lation of Polyakov-type loops for quenched QCD 
is to impose a bias on their value. For exam- 
ple, placing a cut Re ^2 i=x y z Pi < with Pi the 
Polyakov loop in the i-th direction, one selects 
configurations which mimic those of full QCD 
generated with the P boundary condition for sea 
quarks. For the P boundary condition for valence 
quarks we then expect larger values of hadron 
masses. The results plotted in Fig. 3 show that 
the size-dependent shift of masses indeed become 
almost as large as those for full QCD. 

4. Finite-size effect for large volume 

Our results and analyses demonstrate that 
finite-size effects are smaller for quenched QCD. 
This, however, does not mean that a lattice 
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Figure 3. Comparison of quenched hadron masses 
with and without the cut Re ,, , Pi < 0. 

Solid lines are fits to full QCD data from Fig.l. 
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Figure 4. mjf jm p in the chiral limit as a function 
ofl/ia(fm). 



size of order 2 to 3fm often used in quenched 
hadron mass measurements is large enough to 
ignore finite-size errors. This is illustrated in 
Fig. 4 where we plot our Kogut-Susskind results 
at (3 = 6.0 for the ratio mjv/rre p extrapolated 
to the chiral limit as a function of 1/La. The 
GF11 results for the Wilson action obtained with 
a Gaussian smeared sink of a size 0, 1 and 2 at 
(3 = 5.7[8] are also plotted. We observe an al- 
most linear decrease of our data, amounting to 
11% between La « 1.7fm to the largest lattice 
size La « 2.5fm. An 8% decrease of this ratio at 
m q a = 0.01 is enhanced by the m q a — > extrap- 
olation; this indicates an increasingly important 
finite-size effects towards the chiral limit. The 
decrease of the GF11 result extending to a larger 
lattice size La « 3.3fm is somewhat less but still 
significant. We should add that much less size de- 
pendence is found in the GF11 data for the sink 
size 4[8] and also in the previous APE data for 
Wilson action at the same (3 and lattice sizes[9]. 

5. Summary 

We have shown that finite-size effects for small 
lattice sizes are much severer for full QCD. The 
difference originates from breaking of the cen- 
ter Z(3) symmetry due to dynamical sea quarks 
in full QCD, which enhances the amplitude for 
propagation of valence quarks around the lat- 
tice. In quenched QCD the Z(3) symmetry of 



the pure gauge action eliminates such amplitudes, 
and hence leads to a smaller finite-size effect. 

Toward large lattice sizes Liischer's analysis for 
point particles in a finite box leads us to ex- 
pect that finite-size shift of hadron masses be- 
comes exponentially small. Onset of such a be- 
havior is not seen in current quenched data up to 
La « 2.5 — 3.3fm, however. Simulations on sev- 
eral different lattice sizes are necessary to extract 
the reliable value of TOjv/m p at infinite volume. 
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